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Abstract 

We determine the set of the Bloch vectors for N- level systems, generalizing the familiar Bloch 
ball in 2-level systems. An origin of the structural difference from the Bloch ball in 2-level systems 
is clarified. 



*Electronic address: gen@hep.phys.waseda.ac.jp 
Typeset by REVTeX 



1 



1. INTRODUCTION 



The determination of a state on the basis of the actual measurements (experimental data) 
is important both for experimentalists and theoreticians. In classical physics, it is trivial 
because there is a one-to-one correspondence between the state and the actual measurement. 
On the other hand, in quantum mechanics, where a density matrix is used to describe the 




state, it is generally nontrivial to connect them 

It is known that the Bloch vector (coherence vector) gives one 

of the possible descriptions of iV-level quantum state which meet the above requirement, 
because it is defined as a vector whose components are expectation values of some observ- 
ables: For iV-level systems, the number of observables that we need to identify the state is 
since there are N 2 — 1 independent parameters of the density matrix p which is 
Hermitian and subject to Tr p = 1. Actually, if we choose the generators Aj (i — 1, . . . , N 2 — 1) 
of SU(N) (e.g., the Pauli spin operators in 2-level system) for observables of interest, the 
density matrix is determined from their expectation values (Aj) 's: 

N 2 -l 

A = ((X 1 ), . . . , (A^)) - p = h N + ~ ( X A- (!) 

i=l 

(See the next section for details.) Thus, experimentalists have only to measure the values 
of Aj (i = 1, . . . , N 2 — 1) for the identification of a quantum state in the iV-level system. 
(If need be, it is possible to determine the corresponding density matrix by using map 
(JTJ)). The Bloch vector A is defined as A = ((Ai), . . . , (\n 2 -i)) £ M^ 2-1 and this gives 
the desirable description of the states for iV-level systems. However, as concerns the set 
of the Bloch vectors (the Bloch-vector space) not much has been determined so far. (For 
experimentalists, the determination of the set is nothing but to prescribe the range of the 
experimental data observed.) For N- level systems, actually the domain constitutes a subset 
of R^ 2 " 1 but not itself. Notice, however, that not all the vectors A e IR^ 2-1 give density 
matrices by map ([T|). because the positivity of p, one of the indispensable properties of the 
density matrix has yet to be imposed. So far, the complete determination of the domain has 
been done only in 2-level systems: It is a ball with radius 1, known as the Bloch ball. On the 
other hand, for TV-level systems (N > 3), just some properties are known; the Bloch-vector 
space is a proper subset of a ball in R N - 1 ; its 2-dimensional sections are clarified in 3-level 



systems [12] and 4-level systems (2-qubits) 



there appears some asymmetric structure 



[13j; and all of these properties are quite different from those in 2-level systems. The aim 
of the present paper is to determine the Bloch-vector space for arbitrary iV-level systems. 
We also clarify the origin of its structural difference from the familiar Bloch ball in 2-level 
systems. 



2. REVIEW OF THE BLOCH VECTOR FOR 2-LEVEL SYSTEMS AND PREPA- 
RATION FOR ITS GENERALIZATION 



In this section, a brief review of the Bloch vector for 2- 
ration for its generalization to iV-level systems (N > 3) 



evel s yste ms 



ii 



12, 



13j 



and a prepa- 



are given. 



A. Density matrix and Generators of SU(N) 

We begin with the definition and some properties of the density matrix and generators of 
SU (JV). The density-matrix space £ +i i(7Yat) for iV-level systems associated with the Hilbert 
space TCn (— C N ) is given by 

C+, 1 {H N ) = {pe£{H N ) :(i) Trp = l, (ii)/5 = p t , (iii) P< >0(i = l,...,iV)}, (2) 

where C(Hn) (— M(N)) denotes a set of linear operators on H N , pi (i = l,...,N) ith 
eigenvalue of p QjDAs one of the properties of the density matrix, 

(iv) Trp 2 < 1 (3) 

follows from Eq. (J2J) for any p G C +) i(Hn)- Notice that the equality holds if and only if p 
is a pure state (i.e., 3 \ip) e Hn s.t. p = \ip)(ip\). Condition (iv) (in Eq. (jSJ)) plays a role 
of characterizing the Bloch-vector space to be in a ball in - 1 , as will be shown later. In 
the case of iV = 2, conditions (i), (ii) and (iii) (in Eq. (J2J)) are equivalent to (i), (ii) and (iv) 

m 

yp G C(H 2 ), (i) Trp = 1, (ii) pt = p, (iii) Pi > 

4» (i) Trp = 1, (ii) p f = p, (iv) Trp 2 < 1. (4) 

This enables us to characterize the density matrix p for 2-level systems with condition (i), 
(ii) and (iv), instead of (iii). Notice that the sufficient condition (<^=) in Eq. (J3J) does not 
hold for a general case (N > 3). 
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The (orthogonal) generators of SU(N) [3| are a set of operators Aj G C(Hn) (i 
1, . . . , A^ 2 — 1) which satisfy 

(a) A l = \ f , (b) TrA, = 0, (c)TrA,A, 



25, 



,j — (5) 

They are characterized with structure constants fij k (completely antisymmetric tensor) and 
gijk (completely symmetric tensor) of Lie algebra su(N): 

[Xi, Xj) = 2if ijk X k (6a) 

. . 4 

[Xi, Xj) + = — 5ijI N + 2g ijk X k , (6b) 

where [•, •] and [•,•]+ are respectively a commutation and an anticommutation relation. 
(Repeated indices are summed from 1 to A^ 2 — 1.) 

A systematic construction of generators of SU (N) which generalize the Pauli spin oper- 
ators is known 



U CU11SI 



12j; the (orthogonal) generators are given by 
{Ai}£ _1 = {u jk ,v jk ,wi}, 



(7) 



where 



«** = \j )(k\ + \k )(j\, v jk = -i(\j)(k\ - \k)(j\), 
u /■:/■-.■ i ; 

i=i 



z(I ^£(U>0'|-^ + l><Z + l|), 
l<j<k<N, 1 < Z < N- 1, 

with {|m)}^ =1 being some complete orthonormal basis of 7Yat. This gives Pauli spin oper- 
ators (Ai = W12 = o"i, A 2 = t*i2 = 02, A 3 = t&i = <t 3 ) for A^ = 2 with structure constants 



(8a) 
(8b) 
(8c) 



fijk = tijk (Levi-Civita symbol), g ijk = 0, (9) 

and Gell-Mann operators (Ai = m 12 , A 2 = v±2, A3 = Wi, A 4 = m 13 , A 5 = fi 13 , A 6 = w 23 , A 7 = 
V23, Ag = W2) for A^ = 3 with non-vanishing structure constants: 



/123 


= 1, 










^458 


= ^678 


= Vs/2, 








/l47 


= ^246 


— f'257 — 


^345 


— —fl56 — ~ /367 


= 1/2, 


S'ns 


= #228 


— #338 = 


~9at 


* = V3/3, 




#448 


= #558 


— #668 — 


9778 


= -V3/6, 




#146 


= #157 


— 9256 — 


9344 


= #355 = —9247 — 


~#366 



(10) 



Any generators {A'j}^ -1 which satisfy Eqs. © with structure constants f[,- k , g'^ k are con- 
nected with this specific generators (J7J) by some orthogonal matrix V G 0(N 2 — 1) with 
X'i = VijXj and 

fijk = VilVjmVknflmn, 9ijk = Vift ^wVhnQlmn- (H) 

Since the Levi-Civita tensor has a rotational invariance, i.e., VuVj m Vk n ^ijk — det Ve^ = 
ie^fc, the structure constants of SU(2) are limited to = ±£ijk and = 0, while in 
the case of N > 3 there is no rotational invariance of structure constants. 

Finally we notice that the generators Aj's of SU (N) with an identity operator In G C(1~Ln) 
form a complete orthogonal basis of £(7ijv), in the sense of the Hilbert-Schmidt product. 



B. The Bloch vector for 2-level systems 

We review the familiar Bloch vector for 2-level systems (N = 2). From the properties (0) 
for iV = 2, it is easy to show that the following statements (I2) and (II2) hold: 

(I 2 ) Any operator p G £(0.2) with conditions (i) and (ii) can be uniquely characterized 
by a 3-dimensional real vector A = (Ai, A2, A3) G M 3 as 

P = \% + (12) 

(II 2 ) By imposing condition (iv) on the above operator p, the length of A is restricted to 
be less than or equal to 1: 

|a| = v^a7<i- (is) 

From these statements and relation (@J), any density matrix in 2-level systems turns out to 
be characterized uniquely by a 3-dimensional real vector where the length satisfies Eq. (|13j). 
Therefore, if we define the Bloch-vector space 5(R 3 ) as a ball with radius 1: 

B(R 3 ) = { A = (A 1; A 2 , A 3 ) G R 3 : |A| < 1 }, (14) 

its element gives an equivalent description of the density matrix with the following bijection 
(one-to-one and onto) map from £?(IR 3 ) to £ +1 (7i 2 ): 

A -> p = h 2 + -\ifri. (15) 

Z?(R 3 ) is called the Bloch ball, its surface the Bloch sphere and its element the Bloch vector. 
Since the equality in Eq. (fTB^) (i.e., |A| = 1) originates from the one in condition (iv), the 
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surface of the ball (the Bloch sphere) corresponds to the set of pure states and its inside to 
mixed states. The inverse map of (JT5j) is 

p^\i = Trp&i (i = 1,2,3). (16) 

Since we can consider <3Vs to be some observables, the physical meaning of a component of 
the Bloch vector is an expectation value of af X = {(o~i), (02), (03)). An important point is 
that a state can be characterized by expectation values of <7j's which are directly observed 
in experiments. Furthermore, the Bloch vectors allow us to grasp characteristics of states 
from a completely geometrical stand point: The components of the Bloch vector themselves 
give expectation values of <7j's. In addition, for any observable 6 = aS.2 + bi&i (a, bi G R), 
its expectation value is a plus the new z' -component of A multiplied by |6| if we consider 
the direction of b as the new z'-axis (See Fig. [TJ; the probabilities to observe eigenvalues 
a+ \b\ and a — \ b\ of 6 correspond to the normalized rates of P — S- and P — S + respectively, 
where P is the projective point from A to a direction of b and S+ (£-) the points on the 
surface of a positive (negative) direction of b. The specific states such as the Gibbs state 
are also characterized in £?(R 3 ); the Gibbs state pp = exp(—(3H) with Hamiltonian H = 
and temperature 1/(3 corresponds to a line segment between origin O (infinite temperature 
: 1/(3 = 00) and S- (zero temperature : 1//3 = 0). The decohered state after observation 
of 6 corresponds to point P. The time evolution of a state can also be captured as an 
orbit in the Bloch ball and it gives a clear visualization for the Larmor precession process, 
thermalization process or decoherence process and other time evolutions. 

C. Preparation for the generalization of the Bloch vector 

The generalization of the Bloch vector to iV-level systems (N > 3) can be done similarly 
to the case of 2-level systems to some extent. If we choose the generators of SU(N) as the 
observables of interest, then the same statements (Ijv) and (Hat) for iV-level systems hold 
from properties (jHJ): 

(Ijv) Any operator p G £(7Yjv) with conditions (i) and (ii) can be uniquely characterized 
by a (iV 2 — l)-dimensional real vector A = (Ai, . . . , Xn^-i) G M^ 2-1 as 

9 = iV ijV + 2 XiXi ' 
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Sphere : Pure state 
Inside : Mixed state 



FIG. 1: Bloch ball in 2-level systems 

(Hat) By imposing condition (iv) on the above operator p, the length of A is restricted to 
be less than or equal to \/2(N — 1)/N : 



V^<^P>. (18) 



N 

However this does not complete the generalization of the Bloch vector, since properties (i), 
(ii) and (iv) are only necessary conditions for (i), (ii) and (iii) (i.e., to be the density matrix), 
but not sufficient for N > 3. This implies that the Bloch-vector space ^(IR^ 2-1 ) (JV > 3) is 
a proper subset of a ball ([18]). Actually, it has been shown that an angle between any 
two Bloch vectors X a , Xp satisfies 

l>cosZ(a,/3)>--^-p (19) 

if the corresponding density matrices p a and pp are pure, i.e., p 2 a = p a , pi = pp, from the 
following property 

1 > TrpaPp > 0. (20) 

This shows ^(R^ 2-1 ) occupies a relatively small part of the ball and implies an asymmetric 
structure for the case of N > 3. Furthermore, 2-dimensional sections of BfR^ 2 " 1 ) for N = 3 



12j and N = A (2-qubits) [13J are examined in detail, which also support these features. 
However, the complete determination of the Bloch-vector space for an arbitrary dimension 
has not been done. In the next section, we will completely specify it for any iV-level systems. 
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3. THE BLOCH VECTOR FOR iV-LEVEL SYSTEMS 



It is clear that a characterization of the Bloch vector depends on how to integrate condi- 
tion (hi) (instead of (iv)) into the form of Eq. (J17|) . For this purpose, the following lemma 
will be useful. 

Lemma 1 Consider an algebraic equation of degree N > 1: 

N N 

^(-l)X^ = l[(x - Xi )=Q (a = 1), (21) 

3=0 i=l 

which has only real roots Xi G R (i = 1,...,N). The necessary and sufficient condition 
that all the roots x% 's to be positive semi-definite is that all the coefficients a; 's are positive 
s emi- definite: 

Xi > (i = 1, . . . , N) at > (i = 1, . . . , N). (22) 

This lemma can be considered as a corollary of the famous Descartes' theorem (rule's of 
signs) ^?|. (For the reader's convenience, we give a direct proof of Lemma ^ in Appendix 
EJ). From this lemma we obtain 

Theorem 1 Let cij(A) 's be coefficients of the characteristic polynomial det(xl/v — p) where 
p is an operator of the form J77| ) and define 

B(R n2 ^ 1 ) = {A G M^ 2 - 1 : ai (X) > (i = 1, . . . , N)}. (23) 

Then a map: 

A G BiR^- 1 ) ^p = h N+ e C{H N ) (24) 
is a bisection from ^(IR^ 2-1 ) to the density-matrix space C+^TCn)- 

Proof of Theorem [T] 

The operator p in Eq. (|24|) clearly satisfies conditions (i) and (ii) (Hermiticity) from (Ijv). 
By taking account of a fact that an Hermitian operator has a real spectrum, the eigenvalues 
Pi's of p turn out to be all positive semi-definite, i.e., condition (iii) holds from the sufficient 
condition of Lemma E] and the conditions Oj(A) > (i — 1,...,N) in Eq. (|23|) . There- 
fore, the operator is a density matrix and Eq. (J24|) is a map from B(R N 1 ) to C +A (H N ). 
The injectivity (one-to-one property) comes from linear independence of Ijv and Aj's; the 
surjectivity (onto property) comes from the necessary condition of Lemma H 
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Q.E.D 

Theorem^states that the Bloch-vector space for iV-level systems is nothing but ^(IR^ 2 ^ 1 ) 
in Eq. (|23|) with the bijection map (|24j) . which gives the corresponding density matrix. 
Components of the Bloch vector can be considered as expectation values of Aj's 

Aj = Tr p\i (i = l,...,N 2 -1), (25) 

which give the inverse map of (|24p. This means that the state of iV-level systems can be 
completely characterized with (N 2 — l)-expectation values of observables Aj's. 

Notice that the coefficients cij(A) (i = 1, . . . , N) can be written down explicitly if we take a 
matrix representation of the operator Ijv/iV+ AjAj/2. Therefore, the expression of ^(M^ 2 " 1 ) 
in Eq. (J23j) is a practical description of the Bloch-vector space. However, it will be further 
convenient and instructive to express the conditions cij(A) > (i = 1, . . . , N) without resort 
to particular matrix representation. Such an expression clarifies the relation to the previous 
discussion in Sec. 12 CI and the difference between N = 2 and N > 3 cases. To obtain it, we 
use the famous Newton's formulas [12( which connect coefficients (i — 1, . . . , N) (ao = 1) 
and the sums of the powers of roots Xj (i = 1, . . . , iV) of the algebraic equation of iV degrees 
(|21|): Newton's formulas reads 

k 

ka k = ^(-ly^C^a^g (1 < k < N), (26a) 

9=1 

where C N:g = x\. Explicitly, 

a i — C\, 

a 2 = (Cl - C a )/2, 

a 3 = (C 1 3 -3C 1 C 2 + 2C 3 )/3!, 

a 4 = {C{ - 6C 2 C 2 + 8dC 3 + 3C 2 - 6C 4 )/4!, 

a 5 = ■ • • , (26b) 

where Cjq >q is simply denoted as C q here. Applying this to the characteristic polynomial of 
p of the form (fTTj) . one obtains 

l!ai = 1, (27a) 
2!a 2 = 1-Trp 2 , (27b) 
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and 



3!a 3 = l-3Trp 2 + 2Trp 3 , 

4!a 4 = l-6Trp 2 + 8Trp 3 + 3(Trp 2 ) 2 -6Trp 4 , 

5!a 5 = • • • , 



(27c) 



where use has been made of Cjv, g = Trp 9 and CV,i = Tr p = 1. Equation (j27aj) means that 
the condition ai(A) > trivially holds, while Eq. (|27bjl tells us that the condition 02(A) > 
is equivalent to condition (iv) in Eq. 0. The latter has made the Bloch- vector space be in 
a ball (JTHJ). (See statement (IIjv) in Sec. 12 CI ) Consequently we understand that while in 
2-level systems the Bloch- vector space is exactly a ball itself because the coefficients exist up 
to 02, in AT- level systems (N > 3) there are additional conditions 03 > 0, 04 > 0, ■ • ■ , > 0, 
which restrict the Bloch-vector space to be a proper subset of a ball. 

One can further obtain the concrete expressions of coefficients Oj(A)'s in Eq. (J2HJ) in terms 
of the structure constants: 



l!oi 
2!o 2 

3!o 3 

4!a 4 

5!a 5 



,N 



N T 1 " 
(JV-l)(JV-2) 3(iV-2) 



AT2 



IN 



A| 2 + — S'ijfcAjAjAfe 



(iV-l)(iV-2)(JV-3) 3(JV-2)(JV-3), 2 , 3{N 



N 3 



N 2 



\X\ + 



AN 



2(N-2) 3 

H 9ijk\^j^k ~ -^9ijk9klmK*jMK, 



(28) 



where the completely antisymmetric property of has been taken into account to evaluate 
the terms in which appears. See Appendix |B] for the detailed calculations. (Notice that 
dj's in Eqs. have meaning only for % < N.) It deserves to say that since the structure 
constants of SU(N) (N > 3) have no rotational invariance, neither do these conditions. 
Thus the Bloch-vector space has an asymmetric structure in M^ 2-1 for N > 3. 

In the following, we illustrate the Bloch-vector space for 3-level systems, the simplest but 
non-trivial case. For 3-level systems, Eqs. ()28)1 read: a-y = 1, 2!a 2 = (I — tMAI 2 ), and 



3!a 3 



2 
L9 



^l^l 2 ~^9ijk\^j^k 



(29) 
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2/V3 










/ /2/V3 




V \ 


/r / 
























\ i-\ 
\v ^\ 








A3 = 


3 ,2 2 

2 A ' 


















-yvs 






-2/V3 












-2/V3 


A, 


£ {A,,A 2 , 


A3}. Aj = Ag 




A, e {A 4 ,A 


;A 6 




= ±A 3 


A, 6 {A 4 ,A 5 ,A 


6, A7}, \ 3 


= A 8 


Other combinations 



FIG. 2: 2-dimensional sections of the Bloch-vector space for 3-level systems: Large circles are 
sections of the ball |A| < 2/\/3- Grey parts are the domain of the Bloch vector. 

and the Bloch-vector space is a ball in R 8 with radius 2/v3, subject to an additional 
condition 

36 - 9|A| 2 + 9g ijk \i\ 3 \ k > 0. (30) 

By using the explicit values of the structure constants for the specific generators (JJJ), the 
condition (JBTjj) (a 3 > 0) reads 



(31) 



S + 18|A| 2 - 27A 3 (A 2 + A 2 - A 2 - A 2 .) + 6\/3A= 
9V3 {2(A 2 + A 2 + A 2 ) - (A 2 + A 2 + A 2 + A 2 )} 
- 54(AiA 4 A 6 + A1A5A7 + A 2 A 5 A 6 - A 2 A 4 A 7 ) > 



In order to visualize the Bloch-vector space, we use 2-dimensional sections E 2 (z, j) [12, lc 
which are defined as £ 2 (z, j) = {A G B{R 8 ) : A = (0, . . . , 0, A», 0, . . . , 0, Xj, 0, . . . ,0)}. In the 
2-dimensional sections, Eq. (JBTj) are classified into 4 types of conditions: 

Type I: Where Aj G {Ai, A 2 , A3} and Xj = As, then 



A 8 < 1/VS, A 8 > ±V3Xi - 2/V3. 

Type II: Where Aj G {A 4 , A 5 } and Xj = A 3 or Aj G {A 6 , A 7 } and Xj = — A 3 , then 

A,<-, A,->-A 2 --. 
J ~ 3 3 ~ 2 4 3 



Type III: Where Aj G {A 4 , A 5 , A 6 , A 7 } and Aj = A 8 , then 

2 

< 1. 



16 



x/3 X 



(32) 



(33) 



(34) 
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Type IV: Otherwise, 

A? + A?<(0 2 . (35) 

Combined with the ball condition |A| < 2/ a/3, we illustrate the 4 types of sections of the 
Bloch- vector space in Fig. |21 It is clear that the Bloch-vector space is a proper subset of the 
ball. In addition, one sees quite asymmetric structures for Types I, II and III, which stem 
from the absence of the rotational invariance of the generators in the condition (J3U)) . 



4. CONCLUSION AND DISCUSSION 



We have characterized the Bloch-vector space for arbitrary iV- level systems as that pre- 
scribed by Eq. ()23j) : and their explicit expressions are given in Eqs. ()28j) . The essential 
difference between 2-level and iV-level systems (N > 3) is whether they have conditions up 
to a2 (the ball condition) or more (i.e., 03 > 0, a 4 > 0, • • • > 0). Asymetric structures 
appear in iV-level systems (N > 3) because the structure constants have no rotational 
invariance. 

The classification of the states such as pure or mixed states, separable or entangled 
states 0, 3| can be also arranged by means of the Bloch vector: The pure states correspond 
to the surface of the ball (1181) and mixed states inside. As concerns the separability and 



entanglement 



13, 



1 tup 



22j (positive partial 



e famous Peres' criterion 
transpose) in 2 x 2 or 2 x 3 composite systems j^j]. Using Lemma ^ to check the positivity 
in the partially transposed transformation, we can determine the sets of separable and 
entangled states in the Bloch-vector space |24|. 

Although the Bloch-vector space is completely specified, it does not mean that we have 
established a natural parameterization j2] like in 2-level systems. In 2-level systems, it is 
nothing but the Bloch ball, which can be naturally parameterized by the polar coordinates. 
Considering the complex structures (brought about by the remaining constraints 03 > 0, 04 > 
0, • • •) in iV-level systems (N > 3), the notion of the Bloch vector in higher dimensional 
systems might not be as useful as in 2-level systems. However it is still considered to be 
important to express quantum states in terms of the expectation values of observables, since 
then experimentalists can directly determine the states with the use of their experimental 
results. In the circumstances, the classification of the state in the Bloch-vector space is 
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further meaningful so that they can find whether the state is pure or mixed, separable or 
entangled, etc., with their data. 
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APPENDIX A: PROOF OF LEMMA Q] 

We present a direct proof of Lemma ^ 
Proof of Lemma [T] 

[Necessary condition] F From Vieta's formula which connects roots and coefficients: 

JV 

a i= Yl x h x h ■■■Xn (* = h ■ ■ ■ > N ) ( Al ) 

i<h<h<-<ji 

it is clear that Xi > (i — 1, . . . , N) =>■ > (i — 1, . . . , N). 

[Sufficient condition] FLet Oj > (i = 1, . . . , N) and assume at least one of the roots is 
negative definite. (Without loss of generality, we can put %n < 0). Let us define {a,}^ 1 
as: 

TV N-l 

^2(-l) j a jX N - j = (x- x N ) ^(-lYajxW- 1 *-' (So = 1), (A2) 

3=0 3=0 

then clearly the following relations 

a, = cii + aj_iXjv {i — I; • • • A), (A3) 

in which a at = hold. In the case of % = N in (jA3|) . it follows that a^-i < because 
> 0, a/v = and Xjv < 0. In the case of i = N — 1 in (|A3J) . it follows that a N _ 2 < 
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because ajv-i > 0, ajv-i < and xn < 0. Continuing this deduction successively for 
i = N — 1, iV — 2, ■ ■ • in (jA3|) . it follows that a\ < for i = 2 and finally we obtain 
a\ = &i + a Xiv < 0. However this contradicts one of the assumptions a\ > 0. QED 

APPENDIX B: CALCULATIONS OF a;'S IN EQ. (JUJ) 
From Eqs. (0) and ©, one obtains 

TrAj = 0, Tr AjAj = 25ij, Tr XiXjX k = 2zij k , 

.... 4 
Tr AjAjAfcA/ = —5ij5 k i + 2zij m z m ku 

4 4 

Tr AjAj A/cA;A m "^^ij^klra "J^^lm^ijk 2Zij n Z nko Z i m 

Tr \i\j\k\i\ m \ n = ■ ■ ■ , (Bl) 
in which Zijk = Qijk + ifijk- Then it is easy to calculate Cn,q = Tr p q where p is of the form 

Cn,i = 1) 

Cn,2 = ^(4iV + 2iV 2 |A| 2 ), 

C N , 3 = j^(8N + 12N 2 \\\ 2 + 2N 3 X i X j X k g ijk ), 



mkl 



C NA = — ^ (l6iV + 48iV 2 | A| 2 + 16^^^ + 4iV 3 | A| 4 + 2N%X j X k X l g ijm g, 
Cn,5 = ■ ■ ■ ■ (B2) 

Substituting these expressions for CV^'s in Newton's formulas (jZH)) , we obtain the explicit 
expressions in Eq. (|2*K|) . 
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